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Abstract 


Tne Gaiss-Christoiel formula is proposed for the mmerical 
evaluation of integrais involving tne Chapman distribution, The method 
is applied to calculate tne virtual height of reflection of vertical 
incidence pulses as a function of freauency for botn the ordinary and 
extraordi:.ary rays in the oases of transverse: or longitudinal propaga- 
tion, It is also shown that the Chapman distribution may be anproy'- 
mated very closely by a cosine curve. 

It is indicated now the calculated vaiues of the virtus! 
heignt may be used to analyze experimental data, 
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for Virtual Height - B, Friedman 


3 -Bquation 1 = change "x" to 9X9 in the denominator 
7 -¥irst equation from top = 

Remove exponent *2" for denominator of third term, 
Second equation from top - 


~ 
é 


Add exponent "2" for denominator of second term, 
8 = Equation 20 


Add exponent "2" in denominawr inside square root sign to read 


(2? 2 on(x))? 3 Re Ro 
Equation 28 
Change 5d" to read "dx" for wariable of integration 


\O 
{ 


12= In expression for f(x) for Ly 


Remov e-/ from R(R * R,) whick occurs as factor in denominator, 


16- In last squation 
Change mn to read Ch(x,) 
23~ Interchange labeling of axes. 
26- Interchange headings of second and third colamnzs, 


OS 


an Order to obveia the wortual aolgat of @ given lonosphers layer 
as a function of the frequency, it i# necessary to svaluate certain integrals 
involving a Chapman ionisation distrittion which are obtained from the 
Appleton-Hartree formula, Unfortunately, the required integrale are so com 
piicated that it is impossible to evaluate them analytically. This paper 
presents a method for numerical evaluation of integrals involving a Chapman 
distripation, ; 

The method is used to obtain the virtual height for vertical in- 
cidence as @ function of frequency for three distinct cases, Firet, when 
there is no magnotic field, secomi when the magnetic field is transverse to 
the direction of propagation and third when the magnetic field is longitudinal, 
Both the ordinary and extraordinary rays are considered, The regulte are com= 
pared with thoge of Jasger* who used a series expansion to evaluate these 
integrais only in the case when no magnetic field is preseut. The agreement 
between the two different methods is good, Recently Hacke*® propossd another 
method for evaluating the integrals by approximating the Chapman distritmtion 
ty a double parabola, This procedure permits him to express the integrals in 
terms of arc tangente and Logaritnms. 

Tn this paper the various integrale required for the virtual height 
under different propezation conditions are rediced to a atanderd form in 
Sections 2 and 3, The integrals are then transformed in Section 4 to pexmait 
rapid computation by the Gauss—Christofel quadrature formula, 

It will also be shown that the Chepman distribution can be approxi- 
wated by a cogine curve and that the resulting integrals can be evaluated very 
simply by an adaptation of the Gauss—Christofel quadrature formulae, 

in Section 5 a method developed by Booker and Seaton, and Appleton and 
Beynon for the determination of the heights of maxisum ionization density is 
used to analyse goms experizental Gate, However, instead of assuming the 
ionization was parabolically distributed as Booker and Seaton did, ws ageume 
a Chapman distritmtion. Somes data given hy Marning le enalyged and the resulta 
agree with those obtained by him using a more complicated procedure, 


“Jaeger, J. C., Proc, Phys. Soc. 59, 87 (1947). 
*acke, J. E., Jr,, I.R.B., Proc, 36, 724 (1948) 
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The method is also extended to the case where there is a magneto~ 
tonic split. An analysis of data obtained from Pennsylvania State College 
providee sn illustration of a useful method of checking the validity of the 
Chapman hypothesis, 

In Section 6 tables and graphs of the virtual height as a function 
of frequency are given for several cases, The date for the longitudinal case 
is for 4 few propagation conditions of lccal interest, No data have been in- 
cluded for the casa of transverse propagation, The detailed computational 
work in this report is typical and is intended to be iblustrative rather 
than complete. 


2. Tog Iovegrale for Virteal Eeight 
Sonstasr « pulse composed of waves of a circular frequency, 

w , and wavs uowbay, &, waieh is propagated wertigaily, Ye shall derive 
integral expreseions for the virtuel height of the tousephere layer from 
which this pulse ia reflected, 

. Let #4 be the index of refraction in the ionosphere, Accord- 
ing to the Appleton-Hartree theory,if absorption ta neglected, jp is given by 

& formla;* 


j BA a ae _ = ses ieoaimihialien = 
% + } “¢ . 
be aoe sf bermmaee Hiliedt 5 
at Los) a Lox )* 7 


@,8 * Gharge and mes of electron 
% = gumver of tens ia the atmosphere 


x, * af 4 O,* A, A, o/m 


fre = magnetic permeability ef the atmosnhers 


i 


o/s © fe %, oft 


B, © transverse component of the varth's magnetic field, 


ay * lengitudinal component of the earth's magnetic field, 


Bete that the positive siga gives the ordinary ray, while the negative sign five: 
he extra-crdinany ray, 


ak. 


We assume that N varies with the height ascording to the 
Chapmas aistribation! 


(2) Hew chix) = Box [ia or ~ e* 5/2) 


where Xy ie the maximmm ion density of the layer, sre x iz the heights abore 
the earth in scale units defined ty 


Bb = khelght of mastinun lonization 
h = gotual haight above the earth 
et 
He =*# ®%egle hei ghs, % 
a . nee 


let a) . be the critical freeuency, that is, the highest fre. 


quency which will be reflected hy the donmosphere if the magustic field is 
neglected, We can write . 


2 
(3) xs es pa us = Cx yR™ 
we 2 
ha 
where 
(4) R= wf oe 
We also intredues the aotations 
(5) BR = H/o. BL O/H, 
$0 that 


5- 


We shall new obtain an integral for the time of travel 
of & vertically incident wave. The time, dt, to travel a distance, ds, is 
given by 

at = ds/u 


where u ie the group velocity of the wavas, ‘The wave starts at the erouud 


and travels upward until ths height, kh, at which -* O, where 1% is vefleeted 
downwards. Ths time of travel is then 
k 
f ds/u, 
6 


and if ¢ denotes the velocity of light in vacuo the virtual height is defined 
as 


(7) f / *(g 
7 h' = 6 dsju. =bh + ==. 1) @ 
, en 


The second form is seeded later because ths first integral will mot converge 
when we introduce the Chapman dietritation, 


From the well<kmown formala for the group velocity: 


u® if (4k/dw) 
and the fact that 


(8) v= ol and v= w/k 
#0 that 
ke Wy /e 
we find that 
( ~ dw) . 24ew de) , lb ata "w*) 
9) c/a aw 3 2 
jf dw # a(w) 


be 


Consider, now, the case of transverse propagation so that vy = 0 


° 


Then, for the ordinary ray, equation (1) reduces to 


ee se - Ch(x) 
put Ms tS? eS 
R 
Since Ys does not apoear, this formila is the same as in the case when no 
magnetic field is praesent, 
For the extra-ordinary ray, 
» 12 
: 2 K vy 
(42) Ay, t -k- —t 
be-Zo=¥> 
aon 
= 1. oh(x) = Sanaa 
= — 
Ro n° (a Re = ch(x)) 


Note that in this forsmla ba 4s the only quantity dependent on Wand Ch(x) is 


the only quantity dependent on height. & 


stant of the 


frequency to 


47 38 previously defined, la a con- 
ionosphere distribution; it is the ratio of the gyro-marnetic 


critical frequeacy. 


In order to find the group velocity u by meang of formila (9) we mst 
evaluate 
4 2.2 ata “R°) 
(13) 5 (4 Fwe) = aye RD 
a(w) a(R) 
from formula (4). Wow, for the ordinery ray, from equation (11) 
ee ee 
A 58 R" = Ch(x) 
so that (9) yields 
(14) cfu = 2/ me 
In equation (10) put 
hi = ho +Hx! , ds = H dx 
then 
Xo 
(15) a*, = 3. * j | a 7 ds = = *7 
ot fe) 7 o fe) 
<0 ot 


where x is the value of x for which &= O, that is, 


ra 2 
Ch(x,) R 


Note that x = x, corresponds to the height of reflection while x = 


ponds to the grand, 


-~@ corres- 


For the extra-ordinary rey equation (12) may be written 
Re Ch(x) 


2 2 2 ¢ 
ee me ONCE) eal eee 
my (n°. RE - on(x))* 
and 
ie 
tie B we ‘s RE Ch(x) 
pe ( xt Ry = 1* 2 9 
ako (2°. R= ch(x) ) 
t 
so that 
*o ve R° h(x) j 
(17) x! | ax ) —— | 1 + | -1 flex, 
ade - 2 Kf xt ( x°- Re “ oh(x))* ? 


x, is that value of x for which 


la = 9, 
in which cage equation (12) reduces to 


5 42 5 9 
(RF - on(x)}° = RoR. 


If we solve this sauation, we find 
Ch(x,) = R(RIR 
and since Ch(x) is always positive, it is clear that we mst distinzulsh be- 


tween the cases where R is greater or less than R, 


Suopose, first, that R, > R, then from (12) 


t 


(18) Re + Ch(x) - R° 


(R, > R) 


) 
Re ~ (R* - caix)) © 


~k 


and 
x. = 0 
XE 
whes 
(19) Ch(x) = R(R + RO} (R, >R) 


Note that if R(R + EB) is gréater than unity the extra-ordinary ray will pene- 
trate tha layer. 


Now, consider the came wher R. < R, then from (12) we may 


write 
and - 
fox, =? 
when 
(21) Ch(x) = R(R = Ri). 
If RB, = R, then 
f° xt 2 im 
oss {2® ~ cn(7) 
A xe = O when 
Ch{x) = on? P 
Re-write formila (17) as 
pik 2 
(22) ay Ty + RY T, + x, 
where 


a 
a 

ta 

i~ 

i 
i] 
8 a 
o 

| 
x | 
f 

tw 
—— | 


Ch(x) ax 
(24) tT, = f aA eee 
2 2 Wy 
aa 5 (a a, = ch(=)) 


The case of lemcitudinel propagation ia simpler, In this 


case y, * O aud equations (1) reduse é5 
- . —- » Gh(x 
es a 1. y¥, + > tr 
go that 
+ Ch(x) 
ne we bal Re - Chix) = fey 
and 
2 2 ; 
a(R” yw S) i Ch{x) 
fo & 13 a 
an® aa tale 
Finally 
x 
0 Ch(x) 
(25) x! = [ dx = nf + Bu fad - 1 + x 
L eo =| AL ar = R)°K : 


where = ia that value of x for which M = 0. We rewrite (25) as 


ee 


(26) eh em, GE + ae L aa 
us pi * a= a)'e 2 * 
where 
cD 
(27) we ax | -1] 
x 
(28) bLo= f onda 
¢ I. Au 


Since ehtx) 
Ue el Ye eee eee 
aan 6 + 
R(B = 8) 


ie the same as the value of Ww, in (11) af R® 4s replaced by a(a-e ), the 


% 
values obtained in the calculation of (15) can be re-interpreted to give ths 


values of (27). 


3. Summary of Integrale in Standar’ Form 


The integrals required for computing the virtual height are given 
in Section 2 by equations’ (15), (23), (24), (27) and (28). To facilitate com 


putation these integrals will be transformed to the standard forn 


; x 
una i. f(x) Ch(x) dx 
~ © 


“fi = Chix} Chix, ) 


where x, is the appropriate reflection height and ¥(x) takes ws particular 


form for each condition of propagation. A typicel case for transverse 
propagation of the ordinary fay is odtained as follows: 


From equation (14) and (15) 


x, ar) 2 
7 1 = 1 -¥1 - Ch(x) 
t f =| ax rs dx 
- 00 i, ~ ch(x)/R® 00 1 - on(x)/R® 


x 


P [. I) 
—_ a +4 = ch(x)/2* yf - ch(x)/a® 


fhe integral is in the standard fora if we teks, 


42 ee a 
2° +4 - ch(x)/ae } 


while taking account of equation (16) where (x) is given by Chix, ) = R 


1 
i{x) = aaa Snaee RSLS 5 
F h(x.) Ga +i = Ch(x)Jtn{z,} ) 


2 


We may then sumesrige our results as follows: 


fransverse Propagation 
Ordinary Ray 


x = x + 
ot 6 7, 


2 o2 2 2 
Ay, 2 = OR ~ Ch{x), Ch(x,) « #& 


and t ta ofthe form I with 
1 


w+ fi- crise) 


f(x) = 


oli- 


Extraordinary Ray 


a 
2 2 7 
, ; RO + R°* Ch(x) ; 
1 BR 
= nee te RIA. on(x.) = R(R +R) 
Ag a ae - (a?- Chi (x) IB . 


and qT is of the form I with 


Re Ch(x) 


f(x)= 
RVA +3) hae = Re + Ch(x) fenix,) 


The integral ®, is of the form I with 


"5 
nf + Ch(x) - RE 


eee 
(R® “a = ch(x)) BR - Ro + Ch(x) yon(x,) 


8 
Ro Re = CGh(x) ae 
ee we Aerie R C's Ch(x_) = R(R-R, ) 
xt (n°-ch(x))® ReR® {> *o 


and, ty ia of the form I with 
RY = chix) | 


: R . 
(ae + {OWER foR CoG) fon) iy 


The integral tT, ts of the Form 1 with 


RR - ch(x) = Be 


f(x) = a ee beac ag cpr ocean aaa 


(a "oRe ‘ cn(x)) fa, R + Re- Chix) yon(x,) 


f(x) = 


Longitudinal Propagation 


Soe x tL, =F — Lo 
Efe 1 oat ER 


Ay =t- RaT RY ch(x,) =R(R TR) 


Ly. is of ths form ! with 


f(x) 2 


+ : 
RETEST (1 +0 RETRY ) 
L, ie of the form I with 


f(z) #1 


The upper eign corresponds to the ordinary ray the lower sign corresponds to 
the extraordinary ray. 


=-1}- 


4, Calculation of the isiegrale 
4,1 The Gwise-Christofel Quadrature Foraula 


In this section the Gause-Christofel quadrature formula will 
be applied to evaluate the integrals for virtual height, The virtue of this 
particular method is that by using only un points it gives/ approximation to the 
integrand equivalent to that of a polynomial of degree 2n-1 and the weight 
function can be suitably chosen to take care of the singularity of the inte-~ 
g¥and. The Gauss-Christofel formals isu ae follows: | 

Lat w(x) be a non-negative weight function in the interval 
(a,d) and phe) (n = 0,1,2...) the sesociated polynomials of degree 0,1,... 
respectively, which are orthogonel ani normal ower (2,5) with the waight 
function w(x), that is 


5 as 
Q) [ vy(z) g(x) w(x) a= 5 TO BT 


= & =] pp 
Let git) be an arbitrary polynomial of degree On-l at most, then the Gauss- 
Christofel formula etates that 


» 
(2) fi gix) w(x) dx = 7 A, alx,) 


B 


i= 


where ZrXoo ove H, are the roots of bbe) = 0 and AyoAne *« re ere fixed numbers. 


We reproduce a well-known proof* which will give directly the 
method for computing the h, and the polynomials p,(*). Since g(x) is a poly- 
nomial of degree 2n-1 at most, then by ordinary division : 

(3) g(x) = q(x) p(x) * r(x) 


where o{x) snd r(x) are polynomials of degree n-1 at most, By the orthogonality 
prooerty the p,(%) ia orthogonal to any polynomial of degrees less than n, 80 


5 
{ aix) p,(2) w(x) ax = 0 
J 1 


Also, we may write 


(4) r(x) =a pi(x) +...% 8.) By y(*) 


© Yor another proof end further details on this method c.f. 
Svego, G., Orthogonal Polynomials, Amer, Math. Soc. Colloq. Pub. Vol. MAIL, 
1939, p. 46. 
Frank, P, & von Mises, &,, Differentiai-und Integralgleicmungen der Physik, 
ond EA. Vol. I, p. 394. 


ginee v(x) ie of degree n=] at most, How 


b b b 
(5) | ele) wiz) dx -f v(x) w(x) dx = a, f Pol) w(x) dx, 


Las} 


B 
and P(x) ie @ known constent which way be determined from Py { wix) dx = 1, 
~@ 


In order to find &,. Consider the identity obteined from (3) and (4)3 


Bel 
(6) g(x) = q(x) p (x) + ne app, (x) 
bet z= x, _ x, be the a roots of p, (=) = 0, ther 
Bel ‘ 
W) a(z,) = i apy (x) 
¢ 
which are n liseer equetions for the n unknowns Boe Bye cos Bye Y¥f this is 
scived for as we find . 
a” yA e(z,) 4 
1 
and then if we set ds e ay” Bix) wix) dx we have from (5) 
bo a b a 
(6) J elx) wx) axe ST py atx) fo py(z) wis) ax = SoA, ely). 


fo find the ocefficients A; we may chooses particular polynomials 6,(*) o Let 


p(x) 
ej) = sre yes) 
and substitute in (8), It ie clear that 
(4) =O fork #¢ j and g,(x,) =} 


se that 
b p,,(*) 
(9) ry =] Bia, Manx) w(x) dx, 


The polynemials p(x) may be found by solving a set of linear equations. 
% : 
Let p(s )=5 . » o . then by the orthogonality property 


b 
[ va p,, (=) w(x) dx # 0, $= 0, 1, cee ml 
Py 


and we find that 


4 b 
a » f x 3 w(x) az = 0, Jd * Ol eee R=] 


@15- 


The solution of this set of n eouations will determine the ratios of the nti 
coetlicisuts, Soe Bye oou ee Their actusl values may be determined by using 
tae fact that 

. p(x)” wlz) as © 1, 
i & 
It will be noted from equation (9) that # sonxtaat factor multiplying p, (x) 
will not affect the valuse of the Aye Meither in the proof wor in the applica. 
tions is it required that the pb) be normal, It ie ordimarily couvenient for 
fixing the values of the coefficients bd, to take bei 


4.2 Application to the Integrals of Virtual Height 
Now consider the Llutegral 


_2(x) Chix) dx 


. e — + aan : 
1 = Chiz 7; ch(x,) 
which is the general form of ths integrals daweloped in the orewione section, 
We shall transform it to a form in which the Gause-CGhristofel formale is 


applicable. 
x/2 


Let y*® 6 , then 
2 
Ch(x) = ele ye” f2 dx @ = 2ay/y 
and 2) 
-y 2 
oevaewe f ae x) 
v5 i 


2 


het x" = y® * 2t , then y ay = 2t dt | 


8G 


© 2 f(x) sly“= y 


1/2 ~¥5/? -¢ 
L=272 e'e@ + eat | 
- 


TS 


ho ~ Gaix)/en(x,) 


According to the Gauss—Christofel quadrature formla 
there exist constants cy, t, (k=1,2,...,n) such that 


oo -t* D 
[ e- g(t) dt = 5 gq g(t.) 
° k=} 


if g(t) ie a polynomial of degree 2n-1 or less. We shall use the right side 
of this forsmila as an approximation to the integral on the left in case g(x) 


is not a polynomial. 


Consider the set of polynomials p(t), Pal), coe Of degree 
1,2,... respectively, which are mtually orthogonal on the interval (0, w) 
using the weight function, w(t) = Pa than the constants me © 1 Se icie tl 
are the roots of ptt) = 0* , The firet three are as follows and may be readily 


calculated by the msthod shown in the previous section: 


p(t) = {t¥n <1 ) 
2 ott , en 
p(t) = eb = na. 3 qe 


25 fn 


p(t) = $2 + dpe emt (1 * oy) - lao 


where 
(1~8 )(n-2) 


ayn (n-3) 


The constanta G, are given by the formla (9) of Section 4.1 


20 p(t) 42 
x 
> * ; pitt, e=t,) 7 as 


Using this quadrature formula and remenbdering that 


Ch(x) =y, exp (1 - a) /2 


we find that 
2 > Re a t(x)y¥y - ¥, 


I= Cy. 
Vo yy - One) 
Chix ) t 
where the bracket is to be evaluated at the point t 
“Onfortunately neither the Hermite nor Laguerre polynowiais can bs used so that the 
polynomials ani their zeroes would be available from standard tabulations, For 
the Hermite oolynowials the range of the integral is required to be ~ op to w, 
while for tne Laguerre polynomials the weighting function is e-* , 


~L7a 


For the caga n= 4 


= Lhos6 " : 
cy 0, 4425 ty 0.1887 
co = 0, 3987 t, = 0,642 
ec, = 0,04 t, = 1.7911 
3 50 3 19 


4.3 The Cosine Approximation 


The integrals developed in Section 2 may be evaluated more ouicily, 
althmezh slightly less accaravely, by using a cosine approximation to the Chapman 
distribution. In thia way most of the resulta in the Tables in the Appendix may 
be obtained by simpie slide mile calcolations. 

Let x! be the abscissa of the point of inflection of the 
Chapman distribution. we fit a cosine curve to the Ghapman curve at this point 
and obtain the following: 


(1) Ch(x) = . (1 + cos 7) »O>x >= 2!, 


Wote that Gh(x) is very nearly equal to zero for all values of x such that 
x < @e', We shell aceume it to be identically zero for this set of values, 


The typical integral that must be evaluated is of the form: 


Lg 


Here g(x) 19 related to the function f(x) in Section 3,2 by the relation 
aix) = Ch(x) f(x). If we use the approximation (1), J becomes 


4 
See ail 
j? ate) f ae eae a 
“e 2 cht x, ) - 1 = cos Ser 
Let 


(3) 1 + cos =+ = Chix.) (1 * cos @) 


and substitute in the integral. Then, 


ye - (x) a 
(4) gu ontye Bi) BO 
* J, 2-Onlx J 
= - cos ey 
Ch x, 


This integral in certain cases rednces to an slliptic 
integral but we shall evaluate it by using the general method davsloped 
in Section 4.1. 


fhe weight function w(x) will be taken as anity, wut instead 
of using the ordinary algebraic polynomials we shall uee trigonometric poly- 
nomials, The trigonometric functions orthogonal over the interval O to Ti 
ave cos n W(n = 0,1,2,...), which may in fact be expressed as polynomials of 
degree n in cos (2. Ths seross of the a polynomial are then simply 
ay * = 
In this case the orthogonal polynomials are 


2 i 
Py ¥= cos k ix 


1 
Po” OU 


Ti 


Formula (7) of Section 4.1 represents n equations when j = 1,2,..n. The 
coefficients A may be obtained by the following special device, If the 
sowstions are added together and use is made of the ideatities 


a cos k & n= 0 J = 1,2,000n. 
y=i 
we find B ‘ 
fee N= a g(x,) 
Then by comparison with equation (8) of Section 4.1 the h, can ve eesn to be 
r, = n/n j = 1,2,..68 
The integral is then 
ate > e(x,) 


(5) J # 


yee 


where from (3%) 
L+ cos => = Ch (x) (1 + cos a) 


We may express (5) comple:ely in terms of Ch(x) by using (1), 


For, then on(x, ) 
Ch(x, ) = {1 + cos A) 
and 
2x' fCh(x_) — ae) 
(6) lee eee A ee 


In particular if we take n = 3, then 


= yh } 
Chix, ) 0.933 © (x. 
(7) Ch(x,) = 0,500 Ch(x,) 


Chix.) = 0,0669 h(x.) 


and (6) may be evaluated very quickly, 


: 
: 


5. Applications to Analysis of Experiacntal Data 

fhe tebles and graphs given in this report may be used to analy#e 
experimental date on the varistion of virtual height with frequency, ‘The 
observed ionosphere layer will be completely characterized, on the assump- 
x ths 
maximum ion density of the Layer; an the height In kilometers nbows ground 
-of the maximum ion density; H, the scule height in metersof the layer, 


tion of the Ghupman distribution, if thres qiantities are kmowni ¥W 


Consider, first, the case of no magnetic field, i, may be obtainsd 
from the critical frequency, ts by the well-known formula: 
x of 
(1) y= 


° ete 


H and he may be determined by a method similar to that of Booker and Seaton* 
but using a Chapman distribution of the layer instead of « quadratic distri- 
bution, From Equations 15 of Section 2 thse height im kilometers is given by 


. 
ot F t 
: Br Shy Se hig | ax |< - || 
= 90 ota 
=h + EQR) 
where FY (8) ie a known function of R given by Fig. 2. 


fhe otserved data give h! ae a function of R, Weing Eo. (2) we 
can find ne and H, 
For example take R= 0.8, Q(R) = O end then 


h* (.8 a) =h. 
and since Y(.93) = + 1, (55) wel 
h! (.93 i.) +h = H*h <b! (.556),). 


Another approach is to plot h! vs PR). From Equation (2) we should get a 
gtraizht line with slope H and intereept hoe 


* Booker, H. G. and Seaton, S.L., Intern. Asen. of Terr. Maz. and Elee., 
Washington Assembly, September, 1934, 
See also Appletoa, 5.¥. and Baynon, W.J.G., Proc. Phys, Soc, (London) 52 
518 (19u0), 
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Consider, now, the case of longitudinal propagation. 
ively, thes 


be the critical frequoncies for the ordinary and extra-ordinary rays respect- 


Let a,» él, 
éd Ss 4)e < ATE “8 
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where i is determined by Ho 
2 2 
(4) ae Lay Hi js 
Frem (3) 
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ge Wy wy -as 
(5) on ® ; 4 “oo* == 


and from (4) NH may be determined. From Equation (26) of Section 2, 
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Again, plotting h! against YW (RR) will give h and H, but since h. and 


H may also be cbtained by plotting h! ageinet W o(B,8,), & comparison of 


the values obtained for he end H shovld give an indication of how well the 
Chapman distribeation fite the actial distribntion is the layer, 


7 


22m 


iven by 


method ty analysing some data ¢ 


the 


We shall illuetrate 


the 


express the data is 


262 |27%4 1283 | 300 | 317 | 343 


gacyeles we may 


Ae!) 


3. &9 


bs 


~~ 
w& 


Taking f 


Manning* 


following table 
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if h! is plotted against ((R) 


we observe that all 


in the graph, 
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the first two lie on a straight line, From the graph we 
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Also, the fact that the points lie on a straight 


agrees fairly closely with tha values of 
that the assumption of a Chapman distribution of the 


fhe above value of bh 


which was found by Manning. 
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buted according to the Chapman theory, then the above analysis is not valid, 
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When the methods of this section are applied to the analyw#is of 
dats showlng a magneto-ionic split, a series of testa for the validity of 
the Chapman distribution is obtained, Suppose the data for the ordinary 
ray are plotted against YW, (8,8). The first test is whether the result- 
ing curve is a straight line, If it is, the values of he and H are read 


off. Now, sup%ose the date for the extraordinary ray are plotted against 


Walk, B,). The second test is whether this curve is a straight line, If it 


is, the values of hi and H on this curve ars read off, The third test is 
that the values of H and ae for the two curves should be the same, 

These tests were carried out on some experimental data obtained 
at the Ionesphere Propagation Laboratory in Penneylvania State College undsr 
quiet night time conditions. It was asgumed that the propagation was purely 
longitudinal and that no appreciable time delay wes caused by the E-layer. 
The assumption of lowpitudinal propagation at this latitude would iotroduce 
only a slight error®, The rewaits ware very satiefactory and showed that 
the assumptions of the Chapman distribution and the magneto-ionic theory 
wers justified in this case, 

Because this method seama to give a simple way to analyge experi- 
mentel data on the loncaphere, we shall present the analysis of the Penn- 
sylvania State data in detail. The observed h' - f records are shown in Fig, 3. 

The first thing was to find the asymuptotes of ths curve, This was 
done ty considering the ratio Of/ Ah! and sxtrapolating it to sero, We 
find f, = 7.38 a.c., f, 6.54 mc, Then, eolving the equations (3) written 


in terms of frequemcy, 


ie 4 * 4, & 
2 2 
t, ae ts = t t 


we find 
te = 6,97, ty = 0,848 


* Rydbeck, 0.B.H., A Theoretical Survey of the Bowsibilities of Determining 
the Distritution of the Free Eisctrons in the Uppsr Atiosphere, 1942, 
Trans, Chalmers Univ. of fech, Ho. 3. 


85 that 
z= 0.848/6.97 = 0.122 and R = £/6.97 


The functions Wi (2B) and W(B,R, ) are either calouleted from the 


formulas on page le & igor found by interpolation from the graphs. There 
are two possible waye of comparing the experimental data with the theoretical. 
One is to plot the values of yi, (GE) unifomly on a horigontal scale and 


above the corresponding values of BR. The experimental deta which give h? in 
terms of f can be used to read h? in terms of R. The valuea of h! are then 
plotted against the corresponding values of R and the result should ve a 

etraignt line. The resalt is shown in Figs, © and D for both the ordinary and the 
extYaordinary rey, 
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The other procedure is to graph h' against R in the let quadrant 
and YR ) againet Kin the 4th quadrant, 


Now for any point 4, on the n* curve find the corresponding 
point ly vertically beneath i: on thes W curve and then draw a horizontal 
line to the vertical axia anda u5° line to the part ay con the negative horiz- 
ontal axis, The point a. with vertical coordinate equal to that of Qs and 
horizontal coordinate equal to that of ly is a point on the desired straight 
line. ‘the values of h) and H may then be read off, 


The following graph (Fig. H) indicates this procedure for the P, 


date in Fig. B, The resulting line in the left quadrant agress with thet 
shown in Fig. G. 
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TABLE [I 
Cosine Approximation to the Chapman Distribution 
Ch(x) = = Fix 
x axp(1-x-e"*) /2 st 2.8 
0 1,0 1,000 
~ 1 29970 -9964 
- .2 -9891L 29858 
= 3 29753 ~ 9682 
alt 29550 @944O 
= 5 09287 09134 
= .6 8949 °8770 
= al 8547 ° 8352 
~ 8 8082 2 1886 
~ 9 0 1558 01378 
wale 6984 6837 
“1.1 6363 © 6270 
~1,2 05712 5685 
=1.3 5041 - 5090 
=1,4 04369 o 4L9 3 
=1.5 o 3712 «3904 
=1.6 » 3085 2 3338 
=1,7 2 2498 2781 
=1,8 21969 0 2263 
=1.9 01506 oL78y 
=2,0 1124 21350 
=2,1 0795 20965 
22,2 20543 20645 
=2.3 29355 9383 
=2,4 .0221 2.0188 
2,5 20130 2 0066 
-2.6 20072 0003 
“267 0037 O01? 
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Table Bo. 3 
a, 
§ + rm 
x ; is * 4 an * a) " x, 
2 - _ Gala) 
a), ali ; | iat ag 
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Py*} aa + &) 
x? 
x Fig. Fig, 4 Fig. ‘ig. 6 
ol = =6. 47 =6. 46 =6,70 ~6. 27 ~5 90 =§,01 
o2@ 4 =3,16 =3,11 =3.25 = 3.26 =3,.20 =3, 01 
3 6-2.17 =2.18 =2, 2h =2, 29 =2,27 =2,10 
ee -1.64 =1,67 -1,70 =1,6f =1,59 
5 1.28 -1,.35 =1, 26 =1, 25 =1618 ~1 0b 
<6 == 85 - 87 = 985 = oB3 = 173 = ofl 
of == 46 = 643 ~ oil - oh = 23 + .3 
12 + ,53 
18 296 
216 1.80 
ef? 2, 22 
ae« +207 + 13 + (18 + 32 + .74 
82 2,20 
68h 2,02 
285 2.52 
855 2075 
a) 1.11 1.32 1.57 2.77 
292 «61,60 2.9 2.92 
093 = 2033 2.78 
29h 92, 82 
9% 3.06 
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1.15 3.69 5B 


1,25 2.79 


NEW YORK UNIVER 
SITY 
INSTITUTE OF MATHEMATICAL SCIENCES 
LIBRARY 
_ ( Washington Place, New York ZN Y. 


OCT 12 ‘96! 
DATE DUE 


GAYLORD PRINTED INU-S-As 


NYU 
EM-17 


Fr¥edman 


Numerical methods for 
evaluation... 


EM-17 


Friedman 
AUTHOR 
Numerical methods for 
TITLE 


evaluationee. 


DATE DUE 


N. Y. U. Institute of 


Mathematical Sciences 
Lo 2 Waverly Pia 
New York 3, N. Y. 
4 Washington Place 


